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ABSTRACT 

The spatial distribution of Kuiper belt objects (KBOs) in 2:1 exterior resonance with Neptune 
constrains that planet's migration history. Numerical simulations demonstrate that fast planetary 
migration generates a larger population of KBOs trailing rather than leading Neptune in orbital 
longitude. This asymmetry corresponds to a greater proportion of objects caught into asymmetric 
resonance such that their resonance angles, 0, librate about values > n (trailing) as opposed to < ir 
(leading). We provide, for the first time, an explanation of this phenomenon, using physical, analytic, 
and semi-analytic arguments. Central to our understanding is how planetary migration shifts the 
equilibrium points of the superposed direct and indirect potentials. Symmetric libration, in which 
<p librates about ^7r, precedes capture into asymmetric resonance. As a particle transitions from 
symmetric to asymmetric libration, if <f> exceeds its value, ip, at the unstable point of asymmetric 
resonance, then the particle is caught into trailing resonance, while if 4> < ip, the particle is caught 
into leading resonance. The probability that the KBO is caught into trailing resonance is determined 
by the fraction of time it spends with (f> > ip while in symmetric libration. This fractional time 
increases with faster migration because migration not only shifts ip to values < 7r, but also shifts the 
stable point of symmetric libration to values > 7r. Smaller eccentricities prior to capture strengthen 
the effect of these shifts. Large capture asymmetries appear for exponential timescales of migration, 
r, shorter than ^10 7 yr. The observed distribution of 2:1 KBOs (2 trailing and 7 leading) excludes 
t < 10 6 yr with 99.65% confidence. 

Subject headings: celestial mechanics — Kuiper belt — comets: general — minor planets, asteroids 



1. INTRODUCTION 

Transfers of energy and angular momentum between 
disks and bodies embedded within them can explain a va- 
riety of observed dynamical architectures. For example, 
the tight orbits of close-in, extra-solar, Jovian-mass plan- 
ets ("hot Jupiters") are thought to be a consequence of 
orbital migration driven by parent disks composed either 
of viscous gas (e.g., Ward 1997) or planetesimals (e.g., 
Murray et al. 1998). Capture of bodies into mean-motion 
resonances is a celebrated signature of migration (see, 
e.g., Peale 1986, and references therein), showcased re- 
cently by the pair of planets orbiting GJ 876 in a 2:1 res- 
onance (Marcy et al. 2001; Lee & Peale 2002). Migration 
finds application in the solar system as well; it promises 
to explain the preponderance of highly eccentric Kuiper 
belt objects (KBOs) trapped in mean-motion resonances 
with Neptune (Malhotra 1995; Chiang et al. 2003ab). 
The outermost planet in our solar system might have mi- 
grated several AUs outward by gravitationally scattering 
planetesimals (Fernandez & Ip 1984; Hahn & Malhotra 
1999), thereby sweeping its exterior resonances across 
the primordial Kuiper belt and filling them with KBOs 
(Malhotra 1995; Chiang & Jordan 2002). 

Notwithstanding the numerous appeals to disk-driven 
migration, it has been difficult to dispel all doubt re- 
garding its relevance. Worries are perhaps most nag- 
gingly persistent in the Kuiper belt, where the abun- 
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dance of resonant KBOs could, in principle, be explained 
by the following stability argument. One imagines that 
the entire belt was once indiscriminately, dynamically 
excited to large eccentricities and inclinations by one or 
more massive perturbers — perhaps Neptune itself, dur- 
ing that planet's high eccentricity phase (Goldreich, Lith- 
wick, & Sari 2004; see also Thommes, Duncan, & Levison 
1999, 2002)— and that only KBOs fortunate enough to 
be kicked into resonances enjoyed the phase protection 
that permits survival for the age of the solar system. 

Chiang & Jordan (2002, hereafter CJ) propose one test 
of the migration hypothesis that helps to break degenera- 
cies of interpretation. They find that in scenarios where 
Neptune's migration is comparatively fast — occurring on 
timescales shorter than 10 7 yr, for their chosen initial 
conditions — more objects are caught into the 2:1 reso- 
nance with mean libration angles greater than tt. This 
asymmetry in the distribution of libration angles mani- 
fests itself as an asymmetry on the sky: at a given epoch, 
more 2:1 resonant objects will appear at longitudes trail- 
ing, rather than leading, Neptune's. This is a predic- 
tion of the migration model that can be tested observa- 
tionally by wide-angle, astrometric surveys such as, e.g., 
Pan-STARRS (Panoramic Survey Telescope & Rapid Re- 
sponse System). Under either the alternative stability 
hypothesis or the hypothesis that Neptune's migration 
occurred over timescales longer than 10 7 yr, we would 
expect equal populations trailing and leading Neptune. 

The capture asymmetry discovered by CJ — witnessed 
also in numerical simulations by Wyatt (2003) in the 
context of debris disks molded by migrating extra-solar 
planets — arose from a purely numerical orbital integra- 
tion. In this paper, we explain the capture asymmetry on 
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physical grounds, using a variety of qualitative, analytic, 
and numerical descriptions. In uncovering the workings 
of the capture asymmetry on a step-by-step, mechanistic 
basis, we gain insight into its applicability to observa- 
tions. Our focus is on the Kuiper belt, but the setting is 
clearly generalizable. 

Terminology in this field can be confusing; we lay down 
some definitions here. The exterior 2:1 resonance is dis- 
tinguished in offering the possibility of multiple stable 
points for the resonance angle, 

= 2A-A N -n7, (1) 

where A and An are the mean longitudes of the test par- 
ticle (KBO) and the planet (Neptune), respectively, and 
w is the longitude of pericenter of the particle. 3 At mod- 
erate eccentricities of the particle, the angle <j> can librate 
(oscillate with bounded amplitude) about a range of val- 
ues. Particles whose libration centers do not equal tt are 
said to be in "asymmetric resonance" or "asymmetric 
libration" (Message 1958; Frangakis 1973; Beauge 1994; 
Malhotra 1996; Winter & Murray 1997; Pan & Sari 2004, 
and references therein). The usual libration of (j> about 
7r is called "symmetric resonance." We explain the phys- 
ical origins of symmetric and asymmetric resonance in 
§2. We refer to the migration-induced preference for ob- 
jects to be caught into asymmetric resonance such that 
they librate about angles greater than tt as "asymmetric 
capture." Note that asymmetric capture does not refer 
merely to capture into asymmetric resonance. We ex- 
plain the mechanics underlying asymmetric capture in 
§3. Connections to observations are made in §4, and a 
summary and outlook for future work are presented in 
§5- 

2. THE ORIGIN OF ASYMMETRIC LIBRATION 

We begin by explaining the origin of asymmetric libra- 
tion in a static 2:1 resonant potential. Our analysis of 
this special case of the circular, restricted, planar 3-body 
problem lays the foundation for understanding the more 
complicated situation in which the planet is migrating. 
We offer two viewpoints: a physical, qualitative descrip- 
tion in terms of impulses imparted to the particle over a 
synodic period (§2.1), and a graphical, quantitative in- 
terpretation using contour plots of constant Hamiltonian 
(§2-2). 

Many of the ideas contained in this section are not 
new; they may be traced in various forms in pioneering 
work on asymmetric resonance by Message (1958), Fran- 
gakis (1973), Beauge (1994), Malhotra (1996), Winter & 
Murray (1997), and Pan & Sari (2004). In particular, 
Frangakis (1973) and Pan & Sari (2004) highlight the 
key role played by the indirect potential. Pan & Sari 
(2004) additionally offer physical explanations of asym- 
metric resonance, concentrating on the large eccentricity 
limit and on how impulsive torques imparted to the par- 
ticle near its periapsis by the central mass and by the 
perturber can balance. We focus instead on the regime of 
small-to-moderate eccentricities. Our methodology owes 
more to that of Peale (1986). We proceed with our ver- 
sion of the facts in the spirit of pedagogy and to establish 

3 We label p:q resonances such that p > q refers to exterior 
resonances (outside the perturber) and p < q refers to interior 
resonances (inside the perturber). This choice runs counter to con- 
vention, but seems commonplace in the Kuiper belt literature. 



the language in which we will describe the dynamics of 
asymmetric resonance when the perturber migrates (§3), 
a subject which seems to have received much less atten- 
tion. 

2.1. Physical Interpretation of the Direct and Indirect 
Potentials 

Asymmetric libration results from the superposition of 
the direct and indirect disturbing potentials. Consider 
the interaction between the Sun, Neptune, and a massless 
KBO in a reference frame centered on the Sun. Denote 
the mass of the Sun by tuq, the mass of the planet by 
mN, the distance between the Sun and the planet by r\, 
and the distance between the Sun and the test particle 
by r. Then the acceleration of the KBO equals 

r = V r (C/ K o P + i?) , (2) 

where 

C/Kep = ^ (3) 



is the Keplerian potential felt by the KBO, and 

Gtiin Gton^n • r 



R — i?dii 



indirect 



(4) 



is the disturbing or perturbation potential due to the 
planet. Here V r is the gradient with respect to r. 
The perturbation potential divides into a direct part, 
Gm-^/\r — tn|, which governs the direct attraction be- 
tween the planet and the test particle, and an indirect 
part, — GmNrN -i"/|i"n| 3 , which accounts for the accelera- 
tion of the reference frame due to Neptune's acceleration 
of the Sun. Note that in an arbitrary accelerating refer- 
ence frame, there is no guarantee that the acceleration of 
a body can be written as the gradient of a potential. We 
are able to do so because the acceleration of our refer- 
ence frame corresponds to the gravitational acceleration 
of the Sun by Neptune. 

We describe the evolution of the test particle's orbit 
inside the exterior 2:1 resonance in terms of its resonance 
angle, 4> = 2A — An — w. The angle <f> is approximately 
the angle of the line drawn through the planet and the 
particle at conjunction, measured from pericenter. This 
geometric interpretation is only approximate because <f> 
is written in terms of mean longitudes rather than true 
longitudes. 

Peale (1986) describes how <f> evolves under the direct 
perturbation, and we summarize his analysis here. From 
conjunction to opposition, the direct perturbation adds 
angular momentum to the KBO through the azimuthal 
acceleration exerted by Neptune. From opposition to the 
next conjunction, the azimuthal acceleration removes an- 
gular momentum. Over a synodic period, the sign of 
the effect of the direct perturbation on cf) can be decided 
by examining interactions near conjunction, when Nep- 
tune and the KBO are closest. Consider conjunctions at 
point A of Figure 1, en route from pericenter to apocen- 
ter (0 < <j> < 7r). The angular momentum removed from 
the KBO before conjunction exceeds the angular momen- 
tum imparted after conjunction for two reasons. First, 
the two orbits are diverging at point A. Second, the dif- 
ference in angular velocities of Neptune and of the KBO 
is smaller before conjunction than after, so that the two 
bodies spend more time close to each other before con- 
junction. The net effect of the conjunction is therefore to 
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L <0 
n >0 



Force on the Sun from Neptune 
in the center-of-mass frame 



What the KBO feels 
in the Sun-centered frame 
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Fig. 1. — Schematic diagram of the effect of the direct disturb- 
ing potential when conjunctions occur at points A and B (not to 
scale). Solid arrows represent direct accelerations of the particle 
by the planet. Open arrows are the azimuthal components of these 
accelerations. The angular momentum and mean motion of the 
test particle are L and n, respectively. Conjunctions at point A 
(0 < <f> < 7r) remove angular momentum from the test particle, 
increase its mean motion, and accelerate (f>. Conjunctions at point 
B result in opposite behavior. After figure 1 of Peale (1986). 



pull the KBO backwards along its orbit, removing its or- 
bital angular momentum. Its semi-major axis decreases, 
and its mean motion increases. 4 As a result, the next 
conjunction occurs later than it would if the planet and 
the test particle were non-interacting; in other words, if 
< 4> < ir, the direct potential increases <f> over its value 
in the absence of the direct perturbation. We can say 
that (j> is accelerated to larger values. Analogous con- 
siderations apply to conjunctions at point B in Figure 1 
(n < <f> < 2ir); Neptune adds angular momentum to the 
orbit of the KBO and <j> decelerates. 

The tendency of the direct perturbation to restore <j> 
towards n is analogous to the effect of a gravitational 
field on the motion of a pendulum, with <j> = n and <j) — 
corresponding to the stable and unstable equilibrium 
points, respectively. 

To understand asymmetric libration, we extend Peale's 
qualitative description to encompass the effect of the in- 
direct perturbation. As is true for the direct acceleration, 
the azimuthal component dominates the evolution of <f>. 
Consider an infinite line connecting the Sun and the test 
particle. We refer to the side of this line toward which 
the test particle is moving as "ahead" of the KBO and 
to the opposite side of the line as "behind." From con- 
junction to opposition, Neptune is ahead of the KBO and 
therefore pulls the Sun ahead of the KBO. As a result, in 

4 This argument neglects changes in the semi-major axis of the 
particle due to the radial component of the direct acceleration; 
these changes are smaller than those brought about by the az- 
imuthal acceleration by of order the eccentricity of the particle. 
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Fig. 2. — The effect of the indirect potential. Open arrows in- 
dicate the directions that bodies are moving in their orbits. Solid 
arrows indicate accelerations. When Neptune accelerates the Sun 
in an inertial reference frame (left), the KBO feels a fictitious accel- 
eration in the opposite direction in the Sun-ccntcrcd frame (right). 
Here Neptune is shown ahead of the KBO; in the Sun-centered 
frame, the indirect torque removes angular momentum from the 
KBO. 



the Sun-centered frame, the KBO feels a fictitious torque 
which removes angular momentum from its orbit (Figure 
2). Likewise, from opposition to the next conjunction, 
the fictitious indirect torque adds angular momentum to 
the KBO's orbit. Integrated over the synodic period, the 
azimuthal component of the indirect acceleration equals 



(T) 



GmN 



sin A6 dt 



(5) 



where A9 = An — / — vj is the angle between the true lon- 
gitude, An, of the planet and the true longitude, / + zu, 
of the particle. For < <f> < tt, Neptune spends more 
time behind the KBO, and (T) is positive; 5 the angular 
momentum of the KBO in the Sun-centered frame in- 
creases, and (f> decelerates (Figure 3). For it < <fi < 2ir, 
4> accelerates. Thus, the indirect potential restores <f> to- 
ward 0. 

Whereas <f> responds to the direct perturbation like a 
pendulum, it responds to the indirect perturbation like 
a metronome. The direct and indirect accelerations op- 
pose one another; they can be balanced, forming stable 
points at values of 4> intermediate between and n, and 
intermediate between n and 2n. This balance underlies 
the phenomenon of asymmetric libration. 

In contrast, for the 3:2 resonance, Neptune spends the 
same amount of time behind and ahead of the KBO re- 
gardless of the value of its (appropriately defined) res- 
onance angle, and the indirect torque averages to zero 
(see §5 for a proof). Consequently, asymmetric libration 
is impossible for the 3:2 resonance. 

2.2. Contours of Constant Hamiltonian 

The long-term evolution of the KBO's orbit when Nep- 
tune's orbit is not varying can be summarized neatly us- 
ing a Hamiltonian formulation. 6 We seek to write down a 

5 The time spent by Neptune forward and backward of the KBO 
does not alone determine the sign of (T); one must examine the 
exact shape of sinAS over a synodic period. However, for the 
2:1 resonance, the shape does not alter the sign of our argument. 
We show in §5 that (T) = -(Gmfj/ajjirci for the 2:1 resonance; 
ci < when < tj> < n, and ci > when 7r < <j> < 2-7T. 

6 We are allowed to use Hamilton's equations with the potential 
V = Ukcp + R because wc can write all of the fictitious forces 
felt in our accelerating frame as gradients of fictitious potentials 
that depend only on the canonical coordinates, not the canonical 
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Neptune ahead 



Neptune behind 




Neptune ahead 



Neptune behind 




Neptune ahead 



Neptune behind 



Fig. 3. — Portions of the synodic period with Neptune leading 
and following the KBO for = n, n < <j> < 2n, and < <j> < tt (top, 
middle, and bottom, respectively). Solid circles indicate starting 
positions and arrow tips indicate ending positions. For example, if 
<j> = 37r/2, Neptune spends more time ahead of rather than behind 
the particle. 



Hamiltonian having no explicit time-dependence so that 
we can plot level diagrams of the Hamiltonian and trace 
the evolution of cj> graphically. We will make these plots 
in two ways, using exact (§2.2.1) and series-expanded 
(§2.2.2) versions of the Hamiltonian. 



2.2.1. Exact Hamiltonian 
We employ the Poincare coordinates, 

A = M + TZ7, 7 = — w, 

and corresponding momenta, 
A 



- y/l-e 2 ), 



(6) 



(7) 



where A is the mean longitude, M is the mean anomaly, 
w is the longitude of pericenter, a is the semi-major axis, 

momenta. As long as this is the case, we can write the Lagrangian 
C = K — V, where K is the kinetic energy and V includes the 
fictitious potentials. 



and e is the eccentricity, all appropriate to the KBO. Fur- 
thermore, \x = Gtoq. These variables are modified from 
the standard set of Poincare variables in order to describe 
the motion of a massless test particle. The Hamiltonian 
reads 

2A 2 



H 



R(t) . 



(8) 



The perturbation R is a function of An(£), which is ex- 
plicitly time-dependent. 

To eliminate the explicit time-dependence in R, we per- 
form a point transformation. Inside the 2:1 resonance, it 
is illuminating to employ the following coordinates: 

</> = 2A- A N (i) - w, <7 = A-A N (i), (9) 

for which the corresponding momenta equal 

r, iv = A-2r. (10) 

We are interested in the evolution of <p (see §2.1), and 
our choice for a will allow us to average over the synodic 
period easily. 
The new time-independent Hamiltonian reads 



H 



H-n N (T + N) 



2(N + 2T) 2 



n N (r + A0- R{<j>,o,T,N) , (11) 



where «n is the mean motion of Neptune. The price 
we pay in eliminating the explicit time-dependence is 
the appearance of an extra term in the Hamiltonian, 
— nj\r(r + N) = — n-^\r x r|. This term renders H equiv- 
alent to the Jacobi constant. This formulation of the 
Hamiltonian will be useful when we allow Neptune to 
migrate in §3. 

We are interested in the behavior of our system over 
timescales longer than the synodic period. We there- 
fore average over the synodic period to obtain our final 
Hamiltonian. The averaged disturbing function equals 



1 

2^ Jo 
Gton 
2tt 



2tt 



R do~ 



271 



r N 



|r N | 3 



da . (12) 



The total averaged Hamiltonian is 



2(N + 2T} 



r -n N (r+AT)-i? avg (0,r,Af) . (13) 



This Hamiltonian does not depend on a. Therefore N is 
a constant of the motion: dH avg /da = —N = 0. 

Figure 4 displays level curves of H avg as well as indi- 
vidual contributions to _ff avg from the direct and indirect 
terms. In Figures 4, 5, and 6, we set ton = 10~ 3 to q . Had 
we used Neptune's true mass (ton = 5.15 x 1O _5 to.0), 
contours of symmetric and asymmetric libration would 
be less pronounced. When the planet is not migrat- 
ing, the test particle's orbit evolves along a contour of 
constant Hamiltonian. The axes are k — ecos0 and 
h = e sin <f>, so that the distance from the center of each 
plot is the eccentricity and the azimuthal angle is <p. The 
semi-major axis and eccentricity of the test particle are 
related by 

q=U JL V (14) 
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-0.2 0.0 0.2 -0.2 0.0 0.2 

k = e cos k = e cos 



Fig. 4. — Contours of constant Hamiltonian for the exact Hamil- 
tonian (-ffavg) and its components. The contour labels give the 
value of — H in the units GrriQ = 1 and <jn = 1. The total Hamilto- 
nian, Keplerian plus coordinate-transformation terms ("base"), di- 
rect perturbation, and indirect perturbation are individually iden- 
tified. For this figure, N = l.l&^f GrriQa^ and = lO~ 3 rri0. 
Note that half of the contours for the indirect potential correspond 
to negative values of -Redirect- 




-0.2 0.0 0.2 -0.2 0.0 0.2 

k = e cos k = e cos 



Fig. 5. — Contours of constant Hamiltonian with R expanded 
to first, second, third, and fourth orders in e (top left, top right, 
bottom left, and bottom right, respectively) for the same values of 
N, on j and mjM used for Figure 4. The contour labels give the value 
of — -Havg in the units GrriQ = 1 and an = 1. The second-order 
expansion compares favorably with the exact Hamiltonian (Figure 
4, top left panel), with small differences in the values of rfr for the 
stable points. In contrast, the third and fourth-order expansions 
poorly represent the exact Hamiltonian. 



for a given N. As e increases, a increases. Thus, the 
distance from the center of each plot also represents the 
semi-major axis of the KBO. 

The "base" Hamiltonian refers to the combination of 
the Kepler and coordinate transformation contributions: 



A* 



2(iV + 2r) : 



-n N (r + iV) 



(15) 



The contours of constant iAbasc do not depend on the an- 
gular position of the planet or of the test particle. The 
Kepler and coordinate-transformation contributions are 
each negative. The absolute value of the Kepler contri- 
bution decreases as a increases, and the absolute value 
of the coordinate-transformation contribution increases 
with a. For a given adiabatic invariant N, |-ffbasc| con- 
tains a local trough. 

The direct and indirect contributions are much smaller 
in magnitude than the base Hamiltonian and thus are 
not important in most regions of h-k space. However, 
in the trough in Tfbase (near exact resonance), the gra- 
dient in iJbase is small enough that the disturbing terms 
significantly alter the shapes of the contours. The di- 
rect term "fills in" the trough near <f> = 0, leading to 
libration about <f> = it. For appropriate values of N, the 
opposing contours from the direct and indirect pertur- 
bations superimpose on the trough in H^ asc to generate 
contours of asymmetric libration: two unstable points 
appear at (f> — ir and <fi = 0, and two stable points ap- 
pear at < <j) < 7T and n < <f) < 2tt. 



2.2.2. Series- Expanded Hamiltonian 

We now approximate the Hamiltonian using a literal 
series expansion of the disturbing function, as presented 
in Murray & Dermott (1999). We must be careful about 
which terms we include. Because the asymmetric islands 
are generated by small changes in the Hamiltonian, their 
shapes are sensitive to artifacts of the expansion. Fig- 
ure 5 portrays the same Hamiltonian featured in Figure 
4 using an expansion to first, second, third, and fourth 
orders in e. The terms that depend on <j> in this expan- 
sion are proportional to cos(j<j>), where j is an integer 
ranging from to the order in e. To first order, only a 
symmetric island can be represented. To second order, 
we obtain a reasonable representation of asymmetric is- 
lands. To third and fourth orders, more islands appear 
that are not present in the exact solution. These arti- 
facts do not appear for all values of N and on- However, 
when the planet migrates (see §3), N remains constant 
while on traverses a range of values. At least some of 
the resulting combinations of and on exhibit artifacts 
if the expansion order is three or four. 

To avoid this problem, we use the second-order expan- 
sion of R in the integrations that follow: 



R, 



avg,2 



G'rr^^ 
a N 



+ f2e 2 + hie cos <j) + / 53 e 2 cos 2<j>) 

1 



2a 



ecos< 



(16) 
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where a = a-^/a and the /j's are functions of Laplace 
coefficients given in Murray & Dermott (1999). We have 
dropped terms that average to zero over the synodic pe- 
riod. The Hamiltonian, expanded to second order and 
averaged over the synodic period, is 



-^avg,2 — -^basc -^avg,2 



(17) 



This second-order expansion generates contours of asym- 
metric libration but no artificial islands in inappropriate 
locations; it provides a good qualitative representation 
of the actual potential. For example, the contours in the 
top left panel of Figure 4 and the top right panel of Fig- 
ure 5 are qualitatively similar, though the stable points 
in the two plots have somewhat different values of <f>. 

Whereas Equation 16 suggests that asymmetric libra- 
tion arises from the second-order term in a power-series 
expansion of the disturbing function, we emphasize that 
such a conclusion is misleading from a physical point of 
view. Asymmetric libration relies on a balance between 
the direct and indirect perturbations; it depends on the 
indirect term not time-averaging to zero. Expanding the 
direct term to infinite order while excluding the indi- 
rect term would not yield asymmetric libration. Beauge 
(1994) illustrates the difficulty of attributing asymmet- 
ric libration to the second-order term by truncating the 
potential for the 3:2 resonance at second order and find- 
ing asymmetric libration where there should be none. 
Equation 16 should be regarded as simply a useful fit- 
ting formula for the true potential. 

3. THE ORIGIN OF ASYMMETRIC CAPTURE 

We now explain the origin of asymmetric capture when 
the planet is migrating. As the planet migrates outward, 
it may capture test particles that begin outside resonance 
into resonance [see, e.g., Murray & Dermott (1999), and 
references therein]. To be caught into asymmetric res- 
onance from an initially circulating orbit, a test parti- 
cle must first be caught into symmetric resonance. This 
may be seen qualitatively using a series of level diagrams 
for the Hamiltonian, and quantitatively by numerically 
integrating the equations of motion (§3.1). We recall 
from Goldreich (1965) that migration induces an offset 
in the stable point of symmetric libration (§3.2). We ex- 
pand upon these two observations to explain the differ- 
ing probabilities of capture into the two islands of asym- 
metric resonance. We elucidate three separate physical 
effects and provide a series of diagnostic numerical inte- 
grations (§3.3). We provide a sample prediction of the 
analytic theory (§3.4), and compare it with numerical 
experiments that solve for the ratio of particles captured 
into libration about the two islands as a function of mi- 
gration timescale and other initial conditions (§3.5). 

3.1. The Transition from Circulation to Asymmetric 
Libration 

Before a test particle can be captured into asymmetric 
libration, it must first be captured into symmetric libra- 
tion. To explain why this is so, we appeal to the contour 
plots discussed in §2.2. Consider the limit in which the 
timescale for migration is long compared to the libration 
period. Then, over timescales short compared to the mi- 
gration time, the particle's orbit evolves approximately 
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Fig. 6. — Sequence of level diagrams of H aV g for fixed N = 
6.2y 'GrriQ AU 1 / 2 and increasing on, marked on the plots in AU. 
Here, m^j = 10~ 3 mQ. A particle whose a is too large to be in 
resonance stays on approximately the same circulating contour 
at a large distance from the origin before resonance capture. As 
Neptune migrates outward, the resonance expands outward to en- 
counter the particle. 



along a contour of constant exact Hamiltonian, H avg . 7 
Over timescales comparable to the migration time, we 
obtain a qualitative sense of the evolution of the parti- 
cle's orbit by considering a progression of contour plots 
corresponding to different values of on- Such a progres- 
sion is displayed in Figure 6, in which panels correspond 
to the same value of the adiabatic invariant, N, but in- 
creasing values of a^. Even when Neptune is migrating, 
the Hamiltonian (averaged over the synodic period) is 
still given by -ff avg ; while tin and i? a vg are now explicit 
functions of time, N is still a constant of the motion. 

In the panels of Figure 6, a larger value for e corre- 
sponds to a larger value for a by conservation of N (Eqn. 
14). As Neptune migrates outward, the libration centers 
migrate outward as well. A particle having too large an 
a to be in resonance circulates in the outer regions of the 
plots, maintaining approximately its same distance from 
the origin prior to resonance encounter. The resonance 
appears to approach such a particle, which will eventu- 
ally find itself on a contour of symmetric libration. Since 
the asymmetric contours are all surrounded by symmet- 
ric contours and the evolution of the potential is smooth, 
a particle in asymmetric resonance must previously have 
occupied a symmetric contour. 

We confirm the evolutionary sequence from circulation 

7 Level diagrams in this paper are drawn for fixed values of 
(In. Level diagrams showing the exact paths of particles in e-<f> 
space when on varies with time cannot be made because then the 
potential is time-dependent. 
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Fig. 7. — Evolution of <j> from circulation to symmetric libration 
to asymmetric libration as the planet migrates outward. In this 
figure, mjM is set equal to Neptune's true mass. The planet migrates 
according to Equation 19 with to = 300 Myr and o = 23.1 AU. 
Time on the plot is measured from t = 311 Myr when the particle 
is characterized by <p = 0, e = 0.01, and a = 37.67 AU (1 AU 
outside of nominal resonance). 



to symmetric libration to asymmetric libration by nu- 
merically integrating the equations of motion for <j). We 



employ the expanded Hamiltonian, H, 



avg,2 ■ 



As discussed 



in §2.2.2, this expansion reproduces the asymmetric is- 
lands in the exact Hamiltonian qualitatively well. Our 
Hamiltonian equals 



#avg, 2 (0,r,A,i) 



2(A + 2r) 2 



-n N (T + JV)- 



GmN 
a N 



/3ie cos (j) + /53 c 2 cos 2<p) — — e cos 



2a 



where e = e(T, N). We prescribe 



= aN <° (£) 



2/3 



a(fi + f 2 e 2 + 
(18) 

(19) 



where on,o and to are constants, so that a = a(r, N, t) 
and «n = n^(t). The /'s are held constant in our com- 
putations. The form of this prescription for Neptune's 
migration will be motivated in the next section. 
The equations of motion, 



dH, 



avg,2 



dH, 



avg,2 



dT 



(20) 



are integrated using the Bulirsch-Stoer algorithm (Press 
et al. 1992). A sample time evolution for <f> is showcased 
in Figure 7, with input parameters listed in the caption. 
The particle initially circulates, then librates symmetri- 
cally, and finally is captured into asymmetric libration. 

3.2. Migration- Induced Offset in Symmetric Libration 

A key ingredient in understanding asymmetric capture 
is the fact that for symmetric libration, planetary migra- 
tion induces an offset, A</>, in the libration center, so that 



Fig. 8. — Migration-induced offset in the stable point of symmet- 
ric libration. Over the first 0.9 Myr shown, the planet migrates; 
over the remainder of the plot, the planet's orbit is fixed. The 
dashed line indicates the offset due to migration, A<fi, calculated 
using Equation 26 at the plotted time of 0.6 Myr; the agreement 
with the numerical integration is excellent. For this integration, we 
drop the term proportional to cos(2<j!>) in J? a vg,2 and set rriN equal 
to the mass of Neptune. The planet migrates according to Equa- 
tion 19, with to = 10 Myr and on,o = 23.1 AU. Time on the plot 
is measured from f = 10.3 Myr when the particle is characterized 
by if> = 0, e = 0.04, and a = 37.67 AU (1 AU outside of nominal 
resonance). 



the stable point lies at <j) = n + A<fi. Goldreich (1965) 
finds that to first order in the eccentricity of the particle, 
symmetric libration is described by the formula 



d 2 (6<b) 



2 dn N 
-7 ocp — 



(21) 



dt 2 dt 
where S(f> — 4> — 7r is small, and 2ir/j is the libration 
period. The term —dn^/dt arises from migration; since 
4> = 2A — An — ru, (f> due to migration equals — An = — "n- 
Equation 21 integrates to 

1 dn^ 



9 sin jt 



7 2 dt 



(22) 



where 9 is the amplitude of libration. Equation 22 de- 
scribes libration of <f> about 7r + A<j>, where the angular 
offset 



7 2 dt 

When Neptune migrates outward, 
magnitude of the offset is 



(23) 

is positive. The 



A<j> = 



4ir T m T 

where 7] is the libration period, T m 
gration timescale, and T a is Neptune's orbital period. 
The libration timescale equals 



(24) 

on/on is the mi- 



TON 



1/2 



(25) 



where C(a) = l/^6a 3 [2/ 3 i(a) - c^ 2 ] w 0.9. Combin- 
ing these expressions, we find that the offset angle is 

3C 2 T to 1 



A(j) : 



4-7T T m m N e 



(26) 
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Fig. 9. — Schematic of the accelerations of <f> for symmetric 
lib-ration. The letters 's' and 'u' mark stable and unstable equi- 
libria, respectively. When Neptune is not migrating, <fr accelerates 
toward tt (left). Neptune's migration accelerates <j> toward larger 
values (center). These contributions combine to yield acceleration 
toward a stable point greater than n, and acceleration away from 
an unstable point less than 2n (right). 



In Equations 25 and 26, e should be evaluated at exact 
resonance. 

We can verify this offset numerically. The Hamilto- 
nian given in Equation 13 with i? avg expanded to first 
order in e reproduces symmetric libration qualitatively 
correctly (§2.2.2). We integrate the equations of motion 
for this Hamiltonian, again forcing Neptune to migrate 
according to Equation 19. 8 By inserting Equation 19 into 
Equation 26, we see that our prescription for migration 
yields an offset, A(f>, that is approximately constant with 
time, modulo the effect of the changing eccentricity of 
the KBO. Figure 8 displays a sample time evolution for 
4> and verifies that this is the case. The center of sym- 
metric libration is shifted above 7r by an amount, Acj), 
that agrees with Equation 26. The input parameters for 
the integration are supplied in the caption to Figure 8. 

Physically, the offset occurs because as on increases, 
the planet's angular speed decreases. Each conjunction 
occurs later than it would in the absence of migration; 
migration accelerates <j> toward larger values. As a re- 
sult, the stable point of (f> increases; see Figure 9 for a 
schematic description. The offset is akin to the shift in 
the equilibrium position of a spring in the presence of a 
constant gravitational field. 

The stable points of asymmetric libration are also 
shifted forward by outward migration. Figure 10 dis- 
plays integrations, using -ff a vg,2, for which the migration 
is turned off after capture into asymmetric resonance. 
The libration centers in the absence of migration relax 
to smaller values of (f>. Figure 11, analogous to Figure 9, 
explains the shifts pictorially. 

3.3. Deciding Between the Two Libration Centers 

Keeping in mind that all libration centers — both for 
symmetric and asymmetric resonance — shift due to Nep- 
tune's migration, we now turn to the question of which 
island of asymmetric libration captures a greater propor- 
tion of particles. We will refer to the asymmetric libra- 
tion center with < 7r as the "leading" center and to the 
center with (f> > ir as the "trailing" center. We choose 
this terminology because in a snapshot of an ensemble 

8 For this integration, we drop only the term proportional to 
cos(2</>) from Equation 16; this permits ease of comparison with 
later integrations in which we restore this term. 



3tt/2 




3tt/2 



0.0 0.2 0.4 0.6 0.8 1.0 1.2 
Time (Myr) 



Fig. 10. — Migration-induced offsets in asymmetric libration. 
Two evolutions of <j> exhibiting asymmetric libration are shown, 
one ending in the trailing island (</> > n, top), and the other end- 
ing in the leading island (<f> < 7T, bottom). The planet migrates 
from time to 0.9 Myr and subsequently its orbit is fixed. Both 
centers of asymmetric libration are shifted to larger values during 
the migration. The migration parameters are the same as those in 
Figure 8 except that here, i? a vg is expanded to second order in e 
and at the plotted time of 0, <f>0 = 6.19 (top) and 6.22 (bottom). 



no migration 

<|> = 71 



effect of migration 

i=7C 
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i = (|> = 

Fig. 11. — Schematic of the accelerations of <j> for asymmetric 
libration. The letters 's' and 'u' mark stable and unstable equi- 
libria, respectively. When Neptune is not migrating, <j> accelerates 
toward the two stable points (left). Neptune's migration acceler- 
ates <j> toward larger values (center) . These effects combine to shift 
the two stable points to larger (f> and the two unstable points to 
smaller (f> (right). The angle ip < n marks the location of one of 
the shifted unstable equilibrium points. Note how the stable and 
unstable points are squeezed together at (f> < n. 
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of resonant particles (see CJ), particles librating about 
the leading center [<p < n) both attain perihelion and 
spend a greater fraction of the synodic period at longi- 
tudes greater than that of Neptune (see Figure 3). As 
a consequence, these particles appear to cluster at longi- 
tudes greater than that of Neptune. Conversely, particles 
librating about the trailing center (cp > tt) cluster at lon- 
gitudes less than that of Neptune. 

We identify three separate factors that determine the 
relative capture probabilities. 

3.3.1. Migration- Induced Offsets Favor Trailing Island 

If, as a particle transitions from symmetric to asym- 
metric libration, <p is greater than some angle ip, then the 
particle is caught into the trailing island, while if <p < ip, 
the particle is caught into the leading island. At first 
glance, one might expect that ip = tt, but in fact ip < tt. 
Figure 11 illustrates how the addition of the acceleration 
of <p due to migration, <p = —tin > 0, shifts the two sta- 
ble equilibria to larger <p and the unstable equilibria to 
smaller <p. The unstable equilibrium point that was lo- 
cated at <p = tt in the absence of migration is shifted by 
migration to a value ip < tt. 

To estimate tp analytically, we perform calculations 
that are analogous to those that yield Equation 21 and 
that incorporate the second-order resonant term in the 
disturbing function. For small 8(p = (p — tt, 

where A(a,e) = 1 - 8e/ 53 (a)/[2/ 3 i(a) - a~ 2 } w 1 - 34e. 
For A > 0, symmetric libration occurs with frequency 
j4 1//2 7 and asymmetric libration is forbidden. The im- 
possibility of asymmetric libration for e < 0.03 (A > 0) 
is not an artifact of our second-order expansion of the 
Hamiltonian. It is true even under the exact Hamil- 
tonian, as can be seen by careful examination of level 
diagrams of the Hamiltonian. 

For A < 0, asymmetric libration becomes possible; the 
unstable equilibrium point, located at <p = tt in the ab- 
sence of migration, is shifted by migration to 

Ip W TT + A _1 A0< TT , (28) 

where A<p is given by Equation 23. The coefficient A~ x w 
-2.8 for e = 0.04. 

The relative likelihoods of capture into the two islands 
depend on the time spent in symmetric resonance librat- 
ing at <p > ip and <p < ip. Migration shifts the unstable 
equilibrium point dividing the two islands backward by 
A~ 1 A<p. Migration also shifts the center of symmetric 
libration forward by A<p. These two effects compound to 
render capture into the trailing island more likely. 

We can verify numerically the value of the transition 
angle, ip. Figure 12 displays two evolutions of <f> that 
start with slightly different initial values for (p. The solid 
horizontal line marks <p — tt, and the dashed line marks 
(p = ip, evaluated using Equation 28 near the time at 
which the particle trajectories appear to diverge. The 
particle that transitions from symmetric to asymmetric 
resonance when <p < ip < tt is caught into the leading 
island, while the particle that makes the transition when 
<p > ip is caught into the trailing island. To the extent 
that particles in symmetric libration spend more time 




o miiLu i 

0.0 0.2 0.4 0.6 0.8 1.0 
Time (Myr) 

Fig. 12. — Migration-induced offset in the unstable point di- 
viding the two islands of asymmetric resonance. Shown are two 
evolutions of <j> with slightly different initial values of <p. The solid 
horizontal line marks <j> = n, and the dotted line marks <f> = ip, eval- 
uated according to Equation 28 at the plotted time of 0.48 Myr. 
If <fi > ip at the moment the particle transitions from symmetric 
to asymmetric resonance, then the particle is caught into the trail- 
ing island (</> > w); if <j> < ip, then the particle is caught into the 
leading island (<p < it). The location of the dividing angle inferred 
from these integrations is well predicted by Equation 28. Note also 
how the symmetric librations near the plotted time of ~0.3 Myr 
are shifted above <j> = tt; this is the same effect documented in 
Figures 8 and 9 and §3.2. For this integration, ra^ is set equal to 
the mass of Neptune. The planet migrates according to Equation 
19, with to = 15 Myr and cin,o = 23.1 AU. Time on the plot is 
measured from t = 15.45 Myr when the particles are characterized 
by <f> = 1.2873 and 1.2876, e = 0.04, and a = 37.67 AU (1 AU 
outside of nominal resonance). 



with ip > ip than with <f> < ip — and note how the symmet- 
ric librations are shifted above <p — tt > ip — we expect 
more particles to be captured into the trailing island. 
Further quantitative analysis is supplied in §3.4. 

3.3.2. Libration Amplitude in Symmetric Resonance 
and Dependence on Initial Eccentricity 

The amplitude of libration exhibited by a particle in 
symmetric resonance depends on the particle's eccentric- 
ity prior to resonance capture, with larger initial eccen- 
tricities producing larger libration amplitudes. Larger 
libration amplitudes in symmetric resonance imply that 
the migration-induced offsets described in §3.2 and §3.3.1 
are less effective at changing the ratio of the times spent 
at <p < ip and <f> > ip. We therefore expect the capture 
probabilities between the two islands to equalize with 
larger initial eccentricities. Figure 13 illustrates the ef- 
fect. Note that larger eccentricities not only produce 
larger amplitudes of libration, but also reduce the mag- 
nitudes of the offsets A(p and tt — ip, as is evident from 
Equations 26 and 28 (the latter for e > 0.03). 

3.3.3. Uneven Libration Rates and Reversal of Capture 

Asymmetry 

We have argued that the ratio of particles captured 
into the trailing and leading islands is determined by 
the ratio of times a particle spends with <f> > ip and 



10 



Murray-Clay & Chiang 



•G- 



■e- 



2tt 



3tt/2 



7T 



tx/2 



2 2 



3tt/2 

tt/2 




i i i i i ■ 




AA/WV 



Initial e=0.01 
1st order 




Initial e=0.01 
2nd order 



■ i i i i i 




Initial e=0.05 
1st order 




}A0 



Initial e=0.05 
2nd order 



0.0 0.1 0.2 0.3 0.4 0.5 0.0 0.1 
Time (Myr) 



0.2 0.3 0.4 0.5 
Time (Myr) 



Fig. 13. — Effect of initial eccentricity on the evolution of (p. Smaller initial eccentricities yield smaller amplitudes of libration and 
greater likelihoods of capture into the trailing island (<f> > 7r). Initial eccentricities equal 0.01 in the left two panels and 0.05 in the right 
two panels. In the top two panels, (f> evolves under the potential in which the term proportional to cos(20) is dropped, yielding symmetric 
libration. The dashed lines indicate the value of n + A0 evaluated using Equation 26 at the plotted time of 0.2 Myr. In the bottom two 
panels, the potential is expanded to second order in e, yielding the possibility of asymmetric libration. The dotted lines indicate the value 
of x[> calculated using Equation 28 at a time of 0.2 Myr. Solid horizontal lines mark <f> = tt. Remaining initial conditions are the same for 
all four panels: the planet has the mass of Neptune and migrates according to Equation 19, with to = 10 Myr and a>j q = 23.1 AU. Time 
on the plots is measured from t = 10.3 Myr when the particles are at <fi = 0.126 and a = 37.67 AU (1 AU outside of nominal resonance). 



with <fi < ip while in symmetric libration. These times 
are determined not only by the range of values spanned 
by <p during symmetric libration, but also by the rate 
at which the particles librate, i.e., as a function of 
(f>. Immediately before the transition from symmetric 
to asymmetric libration, a particle librates more slowly 
when <fi < if) than when <p > ip, opposing the tendency 
to spend more time at <j) > ip as described in §3.3.1. 
Surprisingly, this effect can be large enough to favor the 
leading island for some initial conditions and to reverse 
the usual sense of the capture asymmetry. 

Figure 14 provides a series of snapshots of an ensemble 
of particles beginning in circulation and ending in asym- 
metric resonance. For this integration, we use the more 
realistic migration prescription, 



a N (i) = ON,f - (fflN,f - «N,o)e 



-t/r 



(29) 



where on.o = 23.1 AU is Neptune's initial semi-major 
axis, dN,f = 30.1 AU is Neptune's final semi-major axis, 
and r is a time constant. In circulation (curve 1), parti- 



cles are evenly distributed over an approximate circle, re- 
flecting the rough constancy of <f>. In symmetric libration 
(curve 2), <p changes more slowly near the turning points 
of the libration. This slowing is evident as a clumping of 
particles near extrema of (j>. 

Immediately preceding the transition from symmetric 
to asymmetric resonance (Figure 14, curve 3), the level 
curve of symmetric libration deforms to just surround the 
islands of asymmetric libration into which the particles 
will eventually be caught. The leading island is smaller 
than the trailing island due to migration-induced offsets 
in the stable and unstable equilibria (see Figure 11). De- 
spite the smaller size of the leading island, however, par- 
ticles still tend to clump in its vicinity, reflecting smaller 
values of cp there than near the trailing island. The re- 
sult, after the last contour of symmetric libration finally 
fissions (curve 4) , is a tendency to capture more particles 
into the leading island. 

The greater degree of slowing near the leading island 
may be crudely understood by noting that at (f> < if), the 
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Fig. 14. — Evolution from circulation (curve 1) to symmetric 
libration (curves 2 and 3) to asymmetric libration (curve 4) of a 
collection of 200 particles that initially share the same value of 
the Hamiltonian. Each point represents the phase space position 
of a single particle recorded at 1 of 4 instants. Neptune migrates 
according to Equation 29, with r = 45 Myr. Snapshots are taken at 
times t = 2.75 Myr (curve 1), 4.40 Myr (curve 2), 4.73 Myr (curve 
3), and 5.10 Myr (curve 4). At t = 0, e « 0.01 and a « 37.67 AU 
(1 AU outside of nominal resonance) for each particle. (Since level 
curves of the Hamiltonian are not exact circles even in circulation, 
initial eccentricities and semi-major axes of particles on a level 
curve are slightly different.) Note, in curve 3, the tendency of 
particles to cluster at cj> < it, despite the greater range of values 
accessible at <i> > w. 



unstable and stable points are closely juxtaposed (Fig- 
ure 11). Since both points correspond to <f> = 0, the 
potential in this region of phase space should be rela- 
tively flat. Particles having small <f> near this turning 
point of symmetric libration spend a comparatively long 
time traversing a nearly flat potential. 

3.4. A Sample Quantitative Estimate 

We use our analytic model to estimate the minimum 
migration speed above which capture into the trailing is- 
land is overwhelmingly preferred over capture into the 
leading island. Capture into the leading island is im- 
possible if, at the time of transition from symmetric to 
asymmetric resonance, the libration angle <f) fails to at- 
tain values less than tp, i.e., if the sum of both migration- 
induced offsets is greater than the amplitude of symmet- 
ric libration: 



+ (tt - V) > . 



(30) 



All quantities in this expression should be evaluated at 
the time of capture into asymmetric resonance. Using 
Equations 26 and 28, we re-write the above criterion as 



T m < (1-A- 1 ) 



3C 2 m Q T 



(31) 



47r wn e8 

We may estimate 9 as a function of e in the adiabatic 
limit. In this limit, at the time of transition into asym- 
metric libration, the contour of constant Hamiltonian 
on which the particle lies (the "asymmetric separatrix" ) 
contains an unstable point at <f> = tt. On the asymmetric 
separatrix, the turning points of libration have approxi- 
mately the same value of e as the value at <fi — it- Then 
Equation 18 requires that 

1 



a(/ 3 ie cos tp + .fee cos 20 tp ) - — e cos 



ftp 



a(/3iecos7r + / 53 e 2 cos27r) — — ecos7r (32) 



where tp gives the values of < 



Ptp 



cos 



2a 

at the turning points: 
B(a) 



(33) 

0.058. The am- 



(34) 

the inte- 
0.04 at 



and B(a) = [af 31 - l/(2a)]/[2a/ 53 ] 
plitudc of libration is 

6(e) = \w-<j> tp (e)\ 

For m N /m = 5x 10~ 5 and a = 37.67 AU, 
gration displayed in Figure 14 offers the value e : 
the time of transition from symmetric to asymmetric res- 
onance. Equation 34 then yields 9 w 1.1, consistent with 
the amplitude shown in Figure 14. From these values, we 
calculate that the critical value of T m below which cap- 
ture into the trailing island is overwhelmingly preferred 
is T miCr it <~ 3.8 x 10 7 yr. In the next section on numerical 
integrations, we check the accuracy of this estimate. 

The only remaining parameter for which we have yet 
to supply an analytic expression is the eccentricity at the 
transition from symmetric to asymmetric libration. To 
estimate this, one could work again in the adiabatic limit, 
using the adiabatic invariant, J = §Td<j) (approximately 
proportional to the area enclosed by the path of a particle 
in h-k space in the small e limit). The value of J in 
symmetric libration equals the value of J appropriate to 
the asymmetric separatrix that the particle crosses. One 
could exploit this fact to identify that separatrix and 
its associated value of e. Another improvement would 
be to compute the relative times that (j> spends above 
and below tp when Equation 30 is not satisfied. Such a 
calculation would have to account for the anharmonicity 
of motion in asymmetric libration (see, e.g., Malhotra 
1996), an anharmonicity that is only accentuated when 
the perturber migrates (see §3.3.3). 

We do not pursue these additional computations an- 
alytically, but rather content ourselves with our phys- 
ical, order-of-magnitudc understanding of asymmetric 
capture embodied in Equations 31-34 and proceed with 
more numerical explorations in the next section. We note 
in passing that the resonance capture theory of Henrard 
(1982) and the analytic estimates of capture probabilities 
derived therefrom (Borderics & Goldreich 1984) cannot 
be applied to the problem of asymmetric capture without 
substantial modification. This is because such theories 
are written in the strict adiabatic limit in which the mi- 
gration timescale is infinitely longer than the libration 
timescale. But asymmetric capture relies upon a finite 
migration timescale. In other words, migration-induced 
offsets are absent from the theory of Henrard (1982). 
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Fig. 15. — Dependence on the exponential migration timcscalc, 
t, of the ratio of particles captured into trailing vs. leading asym- 
metric resonance assuming low initial eccentricities (eo as 0.01). 
Horizontal (solid and hatched) bars correspond to a ratio of oo. 
For fast migration, more particles are caught into the trailing is- 
land, and for slow migration, the ratio is near 1. For all timescales 
with a ratio plotted, 100% of particles arc captured from circula- 
tion. For this plot, the mass of the planet is Neptune's true mass, 
and Neptune migrates according to Equation 29, with a^.o = 23.1 
AU and a^,f = 30.1 AU. Each ratio is calculated by following 
200 particles that initially circulate, have the same Hamiltonian 
value, and are spaced evenly in <f>o- The hatched bar corresponds 
to particles having initial semi-major axes of ao ~ 37.67 AU (1 
AU outside of nominal resonance), and the solid bar corresponds 
to ao ~ 43.67 AU (7 AU outside of nominal resonance). Because 
Neptune's migration rate decreases exponentially with time, Nep- 
tune is migrating more slowly when it captures particles that start 
further from resonance. As a result, particles having ao ~ 37.67 
AU are asymmetrically captured over a greater range of migration 
timescales as compared to those having ao ~ 43.67 AU. 



3.5. Population Ratios as a Function of Migration 
Speed and Initial Eccentricity 

We calculate numerically the ratio of captures into the 
trailing and leading islands as a function of migration 
timescale, initial semi-major axis, and initial eccentricity, 
using the migration prescription given by Equation 29. 
For a given r, a set of ri initial values of (a , e , 0o) are 
computed which correspond to the same values of -ff a vg.2 
and N. If particles do not all start on the same Hamil- 
tonian level curve, phase differences accumulate between 
particles and complicate interpretation of the final cap- 
ture ratio. Initial values of ao and eo are chosen such 
that all particles begin in circulation; initial values of </>o 
arc distributed uniformly from to 27r. To achieve suf- 
ficient statistics, we select 77 = 200 or 800 depending on 
the efficiency of capture from circulation. Figure 15 plots 
capture ratios for particles with small initial eccentricity 
(eo f=s 0.01), while Figure 16 (top) supplies results for 
larger initial eccentricity (eo ~ 0.05). In each figure, two 
curves for two different values of ao are delineated, one 
starting approximately 1 AU away from resonance, and 
the other starting 7 AU away. 

All three effects documented in §3.3 manifest them- 
selves in Figures 15 and 16 (top). Shorter migration 
timescales yield greater asymmetries in the capture ra- 
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Fig. 16. — (Top) Dependence on the exponential migration 
timescale, t, of the ratio of particles captured into trailing vs. lead- 
ing asymmetric resonance assuming high initial eccentricities (eo ~ 
0.05). Horizontal bars correspond to a ratio of oo. Fast migration 
yields large asymmetries favoring the trailing island, while slow mi- 
gration yields ratios near 1. The timescale dividing these behaviors 
is shorter for eo ~ 0.05 than for eo ~ 0.01 (Figure 15); see §3.3.2 
for an explanation. (Bottom) Fraction of particles caught from 
circulation. Note that large asymmetric capture ratios correspond 
to timescales shorter than those producing adiabatic capture frac- 
tions. Migration parameters are the same as those for Figure 15. 
Each datum is calculated by following 800 particles that initially 
circulate, have the same Hamiltonian value, and are spaced evenly 
in 0o- For the hatched bar, ao ~ 37.67 AU, and for the solid bar, 
ao ~ 43.67 AU. Neptune is migrating more slowly when it captures 
particles that start further from resonance. 



tios, with preference given to capture into the trailing 
island, because of migration-induced offsets in the sta- 
ble and unstable points of symmetric and asymmetric 
resonance (§3.3.1). For e « 0.01, r < 10 Myr, and 
our chosen ao's, the probability of capture into the trail- 
ing island is overwhelming. Particles that originate from 
smaller semi-major axes are more likely to exhibit asym- 
metric capture because they are caught at earlier times 
when Neptune is assumed to have migrated more quickly. 
At large r, the capture ratio approaches unity from be- 
low, a reflection of the nearly flat potential at <fi < ip 
(§3.3.3). The ratio of trailing-to- leading particles never 
dips below 0.5. 

Our expectation, based on a simple calculation in §3.4, 
that the capture ratio is infinite if T m = a-^ /<zn < 38 Myr 
is largely borne out in Figure 15. The relevant curve is 
the one for a s=s 37.67 AU, for which parameter values 
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best match those assumed for our calculation. The cap- 
ture ratio is infinite for r < 23 Myr, which translates to 
T m = Ta^(t)/[aN.f — a^(t)] < 87 Myr, a critical timescale 
that differs by a factor of 2.3 from our estimate in §3.4. 
The discrepancy could arise from our use of Equation 34, 
which is only correct in the adiabatic limit and for ec- 
centricities that are not small, and our use of Equation 
25, which is only correct in the limit of small libration 
amplitude. 

As initial eccentricities of particles increase, capture 
into symmetric resonance becomes no longer certain. 9 
For eo ~ 0.05, the fraction of captured objects satu- 
rates at 28% for t > 5 Myr and decreases with shorter r 
(Figure 16, bottom). Of those that are captured, strong 
asymmetries in the trailing-to-lcading ratios do not ap- 
pear except at r ~ 1 Myr because the larger ampli- 
tudes of symmetric libration caused by the larger ini- 
tial (free) eccentricities weaken the effects of migration- 
induced shifts (§3.3.2). Moreover, the shifts themselves 
are smaller because of the larger eccentricities. For exam- 
ple, for t = 10 Myr and eo w 0.01, the ratio of trailing- 
to-leading particles can be infinite, while for the same r 
and eo ~ 0.05, the ratio is near unity. The curves for 
e ~ 0.05 are shifted to lower r's compared to those for 
e « 0.01 by a factor of ~8. 

The computations we have just described and numer- 
ical simulations by CJ agree that the characteristic mi- 
gration timescales required to produce strong asymmet- 
ric capture are in the range 10 6 -10 7 yr. Equation 31 
also suggests the same critical timescales. Our calcula- 
tions yield larger asymmetries, however, than the purely 
numerical ones of CJ. For example, Figures 15 and 16 
indicate that at r — 10 6 yr, the capture ratio is infi- 
nite, while CJ report a ratio of 3-to-l. We are unsure as 
to why there are differences, but the computations are 
at heart different: one employs a truncated version of 
the Hamiltonian and time-averages over the synodic pe- 
riod, while the other has no such approximations. While 
the semi-analytic model has served us well in illuminat- 
ing the numerous physical effects underlying what was 
once a purely numerical result, we should probably rely 
on the numerical simulations for actual comparisons be- 
tween theory and observation. 

4. COMPARISON WITH OBSERVATIONS 

Do the observations indicate symmetry or asymmetry 
in the distribution of 2:1 resonant KBOs with respect to 
the Sun-Neptune line? The Deep Ecliptic Survey (Chi- 
ang ct al. 2003ab; Buie et al. 2003; Elliot et al. 2004) 
has developed a classification scheme that identifies res- 
onant KBOs by virtue of their librating angles; an ob- 
ject is deemed resonant if three orbital integrations, each 
lasting 30 Myr and starting with initial conditions lying 
within the 3er confidence surface of possible osculating 
orbits, all yield libration of the same resonance angle. 
By this criterion, 11 2:1 resonant KBOs ("Twotinos") 
are identified, of which 2 librate symmetrically and 9 
librate asymmetrically. Of the asymmetric librators, 2 
librate about 4> > ir (trailing island), and 7 librate about 
(f> < 7r (leading island). The current positions of all 11 

9 Capture is certain if the separatrix dividing circulation and 
symmetric libration forms outside of the circulation trajectory oc- 
cupied by the particle and probabilistic if this separatrix forms 
inside the circulation trajectory (see, e.g., Peale 1986). 



Twotinos are displayed in Figure 17; they are overlaid on 
theoretical snapshots of Twotinos taken from CJ. 

The observed asymmetry — 2-to-7, in favor of the lead- 
ing island — argues against a rapid migration history 
for Neptune. Given the initial conditions assumed by 
C J — most notably initial eccentricities ranging uniformly 
from to 0.05 — the hypothesis that Neptune's migra- 
tion timescale was as short as on/on ~ 10 6 yr predicts 
that the probability that a particle caught into asym- 
metric resonance is in the trailing island is u ~ 0.75 
(CJ), with shorter migration timescales presumably giv- 
ing rise to greater values of u. Given the observed (pos- 
sibly biased — see below) asymmetry of 2-to-7, the prob- 
ability that u > 0.75 (r < 10 6 yr) is 0.04%, where 
we have used the differential probability distribution 
dP/du = (S+l)Sl/[Yl{S-Y)\}u Y (l-u) s - Y (Port 1994, 
pages 264-265). 10 Here S — 9 is the sample size and 
Y = 2 is the number of observed trailing objects. Thus, 
t < 10 6 yr seems unlikely. The probability that u > 0.5 
(t < 10 7 ) is 5.5%. 

4.1. Theoretical Caveats 

One theoretical caveat to the above interpretation is 
that the capture probabilities into the two islands are 
sensitive to the particles' initial eccentricities (see §3.3.2 
and §3.5). The distribution of eccentricities assumed by 
CJ might not be realistic; see, e.g., Chiang et al. (2003ab) 
who present evidence that a fraction of the Kuiper belt 
was pre-heatcd to large eccentricities prior to resonance 
sweeping. If Twotinos today originated from larger ec- 
centricity orbits prior to resonance capture, then migra- 
tion timescales as short as ~10 6 yr would yield capture 
ratios closer to unity. However, objects on initially larger 
eccentricity orbits are less likely to be caught at all by 
the 2:1 resonance. 

If we assume that Twotinos originated from circulat- 
ing orbits having initial eccentricities < 0.05, then slow 
migration, on timescales longer than ~T0 7 yr, is a more 
viable conclusion. Based on our computations, slow mi- 
gration yields 0.33 < u < 0.5. Given the observed asym- 
metry of 2-to-7, the probability that 0.33 < u < 0.5 
(t > 10 7 yr) is 24%. 

A second theoretical caveat is the possibility that ob- 
jects originally caught into one island transition to the 
other over the age of the solar system (Hahn, personal 
communication). We expect such mixing to afflict only 
objects having the largest libration amplitudes. Mixing 
has been reported by CJ on timescales as short as 10 7 yr 
for objects having libration amplitudes of 45° (see their 
figure 8). The Twotinos reported by the Deep Eclip- 
tic Survey do not exhibit mixing on timescales as long 
as 30 Myr, but longer-term integrations have yet to be 
performed. A third caveat is that massive Twotinos — 
objects having a few times Pluto's mass — can deplete the 
two islands differentially (Tiscarcno & Malhotra 2003; 
Tiscareno 2004). Such large objects either have yet to 
be discovered or were ejected early in the history of the 
solar system. 

4.2. Observational Biases 

10 Chapter 23 on random Bernoulli trials contains typos in the 
formulae for dP/du. 
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Fig. 17. — Current positions of 11 known Twotinos discovered by surveys world-wide. Of these 11, 2 are symmetric librators (open 
squares), while 9 are asymmetric librators (solid circles) that we use to constrain Neptune's ancient migration speed. Overlaid are simulated 
snapshots of 2:1 resonant objects taken from CJ; the left panel portrays the outcome for an exponential migration timescale of r = 10 7 yr, 
while the right panel corresponds to r = 10 6 yr. Two asymmetric librators lie in the trailing island, while seven lie in the leading island; a 
de-biased estimate of the population ratio is 3-to-6 and rules out the hypothesis that r < 10 6 yr at 99.65% confidence. 



What about observational biases? Since the 3:2 reso- 
nance admits no asymmetric libration, Plutinos, unlike 
Twotinos, have no choice but to be distributed symmetri- 
cally about the Sun-Neptune line. For this reason, Pluti- 
nos can serve to calibrate detection efficiencies of surveys 
as a function of orbital longitude. Observationally, 35 
Plutinos have been discovered at longitudes leading Nep- 
tune and 25 have been discovered at trailing longitudes 
(see, e.g., figure 14 of CJ). The observed asymmetry is 
in the same sense as for the Twotinos. If we assume 
that Plutinos are equally likely to be trailing or lead- 
ing Neptune, then the probability that the difference in 
the number of Plutinos found trailing and leading Nep- 
tune is > 10 (given a sample size of 60) is 25% according 
to the binomial distribution. This result docs not in- 
dicate strongly whether discovering and astrometrically 
recovering KBOs at longitudes leading Neptune has his- 
torically been easier than at trailing longitudes. 11 If we 
assume a bias and use the Plutinos to de-bias the Twoti- 
nos, then we estimate a de-biased trailing-to-leading ratio 
for Twotinos of 3-to-6 (2.6-to-6.4, rounded to the nearest 
integer). The likelihoods cited above increase to 0.35% 
(r < 10 6 yr), 17% (r < 10 7 yr), and 39% (r > 10 7 yr). A 
small increase in the total sample size of 2:1 KBOs could 
dramatically increase our confidence that fast migration 
did not occur. We look forward to the advent of the Pan- 
STARRS synoptic survey with its anticipated discovery 
of -10 3 Twotinos. 



5. SUMMARY AND EXTENSIONS 

11 Leading longitudes correspond to spring/summer in the 
Northern hemisphere, and trailing longitudes correspond to 
fall/winter. 



Orbital migration of bodies embedded in disks can ex- 
plain observed structures in systems ranging from plan- 
etary rings to extra-solar planets. As a body migrates, 
it can capture others into mean- motion resonances, leav- 
ing a potentially lasting signature of its migration. In 
our solar system, numerous Kuiper belt objects (KBOs) 
exist in mean-motion resonance with Neptune, suggest- 
ing that Neptune may have migrated outwards by sev- 
eral AUs (Malhotra 1995; Chiang et al. 2003ab). Chiang 
& Jordan (2002, CJ) discover by numerical simulation 
that the subset of KBOs inhabiting the 2:1 resonance 
furnishes a unique probe of Neptune's migration history. 
The resonant angle, <j>, of such KBOs can librate about 
values equal to tt (symmetric resonance), greater than 
7r (trailing asymmetric resonance), or less than ir (lead- 
ing asymmetric resonance). These authors find that if 
Neptune's migration occurred on timescales shorter than 
10 7 years (for their assumed initial conditions), more 
2:1 resonant KBOs would have been captured into trail- 
ing resonance than into leading resonance; as a conse- 
quence, more KBOs would be discovered today at lon- 
gitudes trailing Neptune than at longitudes leading it. 
For longer timescales of migration, more equal numbers 
of KBOs would be found leading and trailing Neptune. 
Observational confirmation of more trailing than leading 
KBOs in 2:1 resonance would constitute strong evidence 
in favor of the (fast) migration hypothesis; while cur- 
rent tallies are more consistent with equal populations, 
a definitive census should be possible with the upcoming 
Pan-STARRS survey. 

We have explained the physical origin of this capture 
asymmetry within the context of the circular, restricted, 
three-body problem. Asymmetric libration in the 2:1 
resonance occurs as a result of the superposition of the 
direct and indirect components of the planet's perturbing 
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acceleration (see also Frangakis 1973; Pan & Sari 2004). 
The direct component, produced by the direct gravita- 
tional attraction between Neptune and the KBO, serves 
to accelerate the longitude of conjunction (<p) toward w. 
The indirect component arises from changes in the ac- 
celeration of the KBO by the Sun brought about by the 
latter's reflex motion induced by Neptune. This accel- 
eration, whose effect on the KBO's orbit depends only 
on the difference between the true longitudes of Neptune 
and of the KBO, accelerates <p toward over a synodic 
period. Thus, the indirect and direct perturbations can 
counter-balance each other to produce asymmetric libra- 
tion of 4> about angles intermediate between and it, and 
intermediate between w and 2ir. 

When employing a literal series expansion of the dis- 
turbing potential, the inclusion of an inappropriate num- 
ber of terms can cause spurious asymmetric resonances 
to appear. For the 2:1 resonance, asymmetric libration 
is represented qualitatively well by a literal expansion to 
second order in e but not by expansions to first, third, 
or fourth order. The second-order expansion should be 
regarded as nothing more than a useful fitting formula 
for the true potential. 

Particles caught by an outwardly migrating 2:1 reso- 
nance fall first into symmetric libration before evolving 
to asymmetric libration. At the time of transition from 
symmetric to asymmetric resonance, if the value of <p 
for a particle is larger than a critical angle ip, then the 
particle will be caught into the trailing island of asym- 
metric libration, while if <p < ip, it will be caught into the 
leading island. The relative probability that a particle is 
captured into the trailing rather than the leading island 
is determined by the fraction of its time spent during 
symmetric libration at <p's greater than and less than ip. 

Three factors determine this fractional time spent: 

1. Migration- induced shifts of the stable and unstable 
equilibrium points of the resonant potential. Plan- 
etary migration shifts the stable points for symmet- 
ric and asymmetric resonance to larger values. For 
asymmetric resonance, migration shifts the unsta- 
ble point formerly at it, to a smaller value, ip < n. 
Analytic theory informs us that the angular off- 
sets of the stable point of symmetric libration and 
of the unstable point of asymmetric libration vary 
inversely with migration timescale, T m = an/on- 
The greater the shifts, the more likely it is that 
particles spend more time in symmetric resonance 
at (p > ip and are captured into the trailing island. 
If the sum of both shifts exceeds the amplitude of 
symmetric libration — true if T m < T m , C rit ~ 10 7 
yr — preference for capture into the trailing island 
is overwhelming. 

2. Initial particle eccentricity. The influence of 
migration-induced shifts is greatest when the am- 
plitude of symmetric libration is smallest; the latter 
demands small particle eccentricities in circulation 
prior to capture. Analytic theory also tells us that 
the magnitudes of the shifts themselves grow with 
smaller eccentricities. 

3. Uneven libration rates during symmetric libration. 
The relative time spent is also affected by differ- 
ences in cp over one period of symmetric libration. 



The average libration rate at <p> < ip is lower than 
at (p > tp because the close juxtaposition of two 
equilibrium points at <p < ip flattens the potential 
there. The difficulty with which a particle traverses 
the turning point at <p> < ip counteracts the effect of 
migration-induced shifts and can even reverse the 
usual sense of the capture asymmetry for large T m , 
but not in a way that gives the leading island more 
than a 2-to-l advantage over the trailing island in 
attracting occupants. 

We confirm the results of CJ that the migration 
timescales necessary to capture more particles into the 
trailing than the leading island are between about 10 6 
and 10 yr. For T m shorter than 10 6 yr, the resonance 
becomes increasingly unable to capture particles at all. 
For a realistic migration prescription (e.g., an exponen- 
tial prescription such as Equation 29), Neptune's migra- 
tion rate decreases with time, so particles beginning at 
smaller semi-major axes are captured while Neptune is 
migrating more quickly and are more likely to exhibit 
asymmetric capture. The final asymmetry in the lon- 
gitudes of the particles is sensitive to their initial semi- 
major axes and, more significantly, their initial eccentric- 
ities. 

Asymmetric libration is not unique to the 2:1 resonance 
(e.g., Message 1958; Frangakis 1973; Beauge 1994; Mal- 
hotra 1996; Winter & Murray 1997; Pan & Sari 2004). 
In a p:l exterior resonance, the resonant angle, 

<p = pX — An — (p — l)tn , (35) 

can librate asymmetrically as well. For all other reso- 
nances, including interior resonances (e.g., 1:2), asym- 
metric resonance does not occur. Asymmetric libration 
of the type exhibited by the 2:1 resonance requires a non- 
zero indirect acceleration. We prove below that the in- 
direct acceleration always averages to zero over one syn- 
odic period for a p:q resonance except when q = 1, where 
p > q refers to an exterior resonance and p < q to an 
interior resonance (p ^ q are relatively prime, positive 
integers). The order of the resonance is s = \p — q\, and 
the resonant angle is (p — p\ — qX^ — (p — q)vj- Our proof 
is essentially identical to that of Frangakis (1973), and 
derived independently; we offer our expanded version for 
ease of reference. 

The azimuthal component of the indirect accelera- 
tion of the KBO by Neptune equals — (Gm^/a^) sin A9, 
where A9 = An — / — is the angle between the 
true longitudes of the planet and of the particle. The 
proof for the radial component reads completely analo- 
gously. Over one synodic period, the tangential acceler- 
ation time-integrates to 

(T) = - < ^H±E / sin (A N - / - w )dt 



GmN 



sin(A0) 



-d(A9) . (36) 



*n JO n N - f(A9) - m 
After one full period of the particle, the particle returns 
to the same true longitude while the planet increases 
its true longitude by 2irs/q, increasing A9 by 2irs/q. It 
follows that f(A9) is periodic in A9 with period 2its/q, 
as is any function of /, including 

F{A9)= [n N - f{A9)-ujY 1 , (37) 
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where 21/ = 2-rrs/q. We have neglected the tiny varia- 
tion of vj <C riN over a synodic period. Plugging this 
expression for F into Equation 36, we find 
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Thus, (T) 7^ only when = s = |p — q| for some j. 
This is satisfied only when q = 1 and j = |p — q\ under 
our assumption that p and q are relatively prime. 12 We 
conclude that asymmetric libration, of the kind exhib- 

12 When q is even, one can see that the indirect term integrates 
to zero by appealing to diagrams analogous to Figure 3 and noting 
symmetry over the synodic period. 

13 The only resonance for which our proof does not apply is the 
1:1 resonance. But it is clear that the indirect acceleration docs not 



ited by the 2:1 resonance, can only exist for exterior p:l 
resonances. 13 

Future theoretical work should focus on determining 
the degree of stochasticity in migration driven by scatter- 
ing planetesimals, and on the size spectrum of planetesi- 
mals required for resonance capture to proceed smoothly. 
We have made preliminary estimates which suggest that 
for a Neptune-mass disk in the vicinity of Neptune, the 
disk mass must be concentrated in planetesimals whose 
sizes do not exceed ~100 km for the theory that we have 
described in this paper to apply. Reasons to imagine that 
planetesimal sizes were small can be found in the review 
of planet formation by Goldreich et al. (2004) . 
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time-average to zero in tadpole and horseshoe orbits, which may 
be regarded as asymmetric and symmetric librations, respectively. 
See Frangakis (1973) and Pan & Sari (2004). 
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